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CALCULATED AERODYNAMIC LOADINGS OF M, W, 

AND A WINGS IN INCOMPRESSIBLE FLOW 

• 

By Frariklin W. Diederich and W. Owen Latham 


SUMMARY 


Presented are the results of theoretical Incompressible-flow calcu- 
lations of the spanwise lift distributions, lift-curve slopes, spanwise 
centers of pressure, aerodynamic centers, coefficients of damping in 
roll, rolling-moment coefficients due to aileron deflection, and induced- 
drag coefficients of twenty M, W, and A wings. These results are com- 
pared with similarly calculated results for ordinary swept and unswept 
wings. 


INTRODUCTION 


Interest in M, W, and A wings has recently increased as a result of 
their potentially favorable stability and aeroelastlc characteristics. 

In order to facilitate the aerodynamic and aeroelastlc appraisal of wings 
of this type theoretical incompressible-flow lift distributions have been 
calculated for twenty plan forms of the M, W, and A types and four angle- 
of-attack conditions. These lift distributions have been integrated to 
yield lift-curve slopes, spanwise centers of pressure, aerodynamic 
centers, induced drags, coefficients of damping in roll, and rolling 
moments due to aileron deflection. The results of these calculations 
are presented in this paper and compared with similar results for ordi- 
nary swept and unswept wings. The method by which the spanwise lift dis- 
tributions were calculated and some of its limitations are discuesed. 


SYMBOLS 


a. c. aerodynamic -center position (of wing, unless specified 

otherwise) , fraction of mean aerodynamic chord rearward 
of leading edge of mean aerodynamic chord 
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A 

b 

c 

c 



d 


F 

q 

s 

S 

V 

w 

X 

y 


aspect ratio (b^/s) 
wing span - 

chord, measured parallel to free stream 
section lift coefficient 

average chord 

wing lift coefficient ^L/qS^ 
wing lift-curve slope, per radian 

wing rolling-moment coefficient moi^n t 

\ ^ 

rolling-moment coefficient for linear antisymmetrical angle- _ 
of-attack distribution from 0 at the root £0 1 radian at the 



coefficient of damping in roll 




induced-drag coefficient 1 - - 

distance of quarter-chord point of mean a^rot^n^ic chord 

reaivard of intersection of quarter-choi’d line and plane of * ,, 

symmetry ‘ " 

downwash factor 

dynamic pressure 

semispan of horseshoe vortex = ^ ^ 

wing area = 

free-stream velocity “ 

downwash velocity 

longitudinal ordinate - 

lateral ordinate ^ 
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y* dimensionless lateral ordinate 


y lateral ordinate of spanwlse center of pressure 

y lateral ordinate of centroid of wing area 

oc^ aileron effectiveness parameter 

X taper ratio 

A angle of sweepback at quarter-chord line 

Subscripts : 

B at position of spanwlse discontinuity of angle of sweep 

1 Inner 

MAC pertaining^, to mean aerodynamic chord 

o outer 

r root 

t tip 

V made dimensionless by dividing by s 

5 pertaining to a unit equivalent aileron deflection 

METHOD AND SCOPE OF CALCULATIONS 

« 

Basis of Method 

According to two-dimensional thln-alrfoll theory the lift on a 
section with angle of attack and parabolic camber may be calculated by 
locating a vortex at the quarter-chord point and setting the ratio of 
Its Induced downwash at the three-quarter-chord point to the free-stresm 
velocity equal to the slope of the airfoil at that point. Since Induc- 
tion effects on unswept wings are generally In the nature of Induced 
angles of attack and Induced parabolic cambers, quite accurate lift dis- 
tributions can be calculated on the basis of a bound vortex at the 
quarter-chord line and the slope at the three-quarter-chord line. This 
simple concept may be applied to swept wings as well, except that 
Induced effects at the root are such as to Invalidate the concept In 
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that region. Nonetheless^ methods based on this cpnce'pt have been used 
extensively for calculating lift distributions on swef)t wings, either 
with continuous distributions of vorticity over thp sjSfen (reference l) 
or by means of horseshoe vortices (reference 2 , for instance) . 

In applying the same method to wings of unconven'y.onal plan form,_ 

two difficulties arise. Both of the approaches usfed in references 1 

and 2 imply the assumption that the lift distribution can be approximated 
with sufficient accuracy by the first three or four terns of a Fourier 
series, and both of the approaches consider the dcnmwash at only three 
or four points along the three-quarter-chord line. Neither of these 
approximations may be valid for wings with an unconventional plan form 
or a complicated hngle-of-attack condition; however, the approximations 
can be avoided in part by using horseshoe vortices ; as dn reference 2 
but treating each individually rather than relating them to each other 
by Fourier series. -- ‘ 3 - 

The resulting procedure is straightforward but, nonetheless, more_^ 
time consianing than the methods of references 1 and 2._ It has the advan- 
tage that it can be applied to plan forms and anglfe-o:^ attack distri- 
butions which cannot be treated by those methods. 'Also, the results of 
the procedure are obtained in a fom which lends i^self^to the calcu- 
lation of aeroelastic phenomena. The accuracy of the method is, for 
plan forms to which the methods of references 1 and 2 apply, at least as 
good as that of these methods. However, in applying the method to plan 
forms such as that of figure 1 , the accuracy is affected adversely by 
the fact that the quarter-chord- vortex thiee-quarter-chord-downwash con- 
cept is not valid at the spanwise discontinuities in angle of sweep. 

Since this effect is local and tends to average for the entire wing, the 
resulting accuracy may be adequate for many purposes. As a result of 
these considerations, this method was used in the calculations described 
in this paper. 


Outline of Method of Calculating Lift Distrib.utlons 


The horseshoe vortices are considered to be centered on the quarter- 
chord line (see fig. l), and the downwash is calculated at corresponding 
points on the three-quarter-chord line. In so doing, the charts and 
tables of the downwash due to a horseshoe vortex with iinlt circulation 

at various points relative to the vortex given in reference 3 can be 

used. In these charts and tables the downwash is presented in the form 
of a downwash factor F. The downwash angle at any point due to any 
horseshoe vortex may be determined by means of the equation 


w _ 

V 8 ns 
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where c^c is the loading coefficient at the station at which the horse- 
shoe is centered, s the semispan of the horseshoe vortex, and F the 
factor appropriate to the relative location of the downwash point and 
the center of the vortex. This factor is determined by calculating the 
lateral and longitudinal distances. Ay and Ax, between these points 
and dividing them by s to get dimensionless distances Ay^ and Ax^ 

(in the notation of references 2 and 3)# respectively; the factor corre- 
sponding to these values may then be obtained from the tables of refer- 
ence 3. 

The downwash angle at one point due to all vortices is obtained by 
performing the following summation-. 



(h = 1, 2, . . . n) (2) 


where the subscript h designates the downwash point, j the vortex, 
and n the number of vortices on the total span. This addition may be 
shortened by setting 


^h;3 = FhJ ^h,n-.l-J 





(3) 


^hj ' - ^h,n-i-l-j 


( 4 ) 


where Fj^j is the downwash factor for a vortex on the right wing 

and Fjj^ n+l-J factor for the corresponding vortex on the left wing; 

the addition can be then confined to the right wing alone, so that 



^h = 1, 2, . , . (5) 
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n-l 

:: 

(v)h " ^ ^ 

^h = 1, 2, . . . (6) 


in the symmetric and antisymmetric cases, respectively. (The upper 
limits In the summation apply to an odd number of voirtic^s numbered from” 

the right wing tip inward, with one of them, number on the airplane 

center line; if there Is an even number of vortices, with the legs of 
two vortices coinciding with the center line, the upper Timit In both 
equations (5) and (6) is n/2. ) 

Equations (5) and (6) are, in effect, simultaneous equations 
for ^2^* ^ 2 or n/2 unknown values of c^Cy/cT In terms of the 

corresponding number of known angles of attack along the_.span. The 
values and may be used directly as the coefficients of these _ 

equations disregarding the term c’/Sits; howev^, in that^case the solu- 
tions of the equations must each be multiplied by the factor 8rts/^ in 
order to obtain the corresponding values of Cjcj'c. From these values 
the lift coefficient may be obtained by integrating the .curve c,c/c’ 
over the span. -- ' 

Preliminary calculations were made with several vortex patterns, 
the one of reference 2 shown in figure 1, one with 20 equally spaced 
vortices of width equal to 0. 1 of a semispan, and some wTth more or 
fewer vortices. For instance, lift distributions were calculated by 
means of two of these patterns and by means of the method of reference 1 
for a wing of aspect ratio 4.5, taper ratio 1, with 45° sweepback. These 
lift distributions are shown in figure 2j the dlstriliutlbns and the llft-^ 
curve slopes are in excellent agreement. The possibilit^of using horsed 
shoe vortices with bound parts slanted to follow the quafter-chord line 
more closely was also investigated. ' T 

As a result of these preliminary calculations a rectangular vortex 
pattern with 20 equally spaced vortices reaching from; wing tip to wing 
tip appears to be the optimum arrangement from the point j>f view of 
computing effort commensurate with attainable accuracy. Tn view of the 
inherent limitations of the theory the increase in accuracy obtainable 
by resorting to a greater number of vortices or to the s^e number of 
vortices with skewed bound parts is largely spurious.; The calculations 
of the present paper have been performed by means of this, vortex 
representation. 


m UV' T'lV m ^ \ gr.-^ 
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The factors F were obtained from reference i)-, the values of F 
and F from equations ( 3 ) and (U), and the values of c^c/'c for various 
distributions of angle of attack or w/v by solving equation (5) for 
symmetrical cases and equation ( 6 ) for anti symmetrical cases. The value 
of s is 0.05b/2 and hence "c/a is equal to 4 o/A. 


Calculation of the Lift, Rolling Moment, 

Bending Moment, and Induced Drag 

The calculated lift distributions were integrated numerically using 
the following factors, the lifts being those calculated by solving equa- 
tions ( 5 ) and ( 6 ) and the factors being based on parabolic approximations 
to segments of the distributions and on the assumption that the lift 
distribution goes to zero with infinite slope at the wing tip: 


Lift at 

y* = 


Calculation 

of - 

Total 

lift 

Moment 

(symmetrical 

case) 

Moment 

( ant 1 symmetric al 
case) 

0.05 

0.0958 

0.0047 

0.0044 

.15 

.1042 

.0149 

.0149 

.25 

.1 

.0250 

.0250 

.35 

.1 

.0350 

.0350 

.45 

.1 

.0450 

.0450 

.55 

.1 

.0550 

.0550 

.65 

.1 

.0650 

.0650 

.75 

.1 

.0750 

.0750 

.85 

.0983 

.0832 

.0832 

.95 

.0956 

.0906 

.0906 


The induced drag was calculated in effect by obtaining the value of the 
lift at values of y* equal to O.9808, 0 . 9239 , O.8315, O.707I, 0 . 5556, 
0. 3827, 0. 1951, and 0 from faired curves of the spanwise lift distri- 
bution, calculating the lifting-line downwash frcm these values by means 
of the downwash factors of reference 4 (with m = 15 in the notation of 
reference 4 ), and integrating the products of the lifting- line downwash 
and the local lift by means of the integrating factors of reference 4 . 


Calculation of the Aerodynamic Center 

Neither the method used in this paper for calculating spanwise 
lift distributions nor that of reference 1 furnishes any information 
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concerning the local centers of pressure of these lift distributions and 
hence the wing aerodynamic centers. The frequently made" assumption that 
these methods imply a two-dimensional pressure distribution with a center 
of pressure at the quarter-chord point is not valid.. In this paper the 
wing aerodynamic centers have been calculated on the ba^is of this assump- 
tion partly because the local aerodynamic centers are not known reliably 
and partly because the net effect of local aerodynamic centers ahead of and 
behind the quarter-chord line on the wing aerodynamic c^ter may be small. 
In order to gain some measure of the error introduced in the wing aero- 
dynamic center as a result of this approximation winjg aa'todynamic centers 
have also been calculated for some plan forms on the ba^s of assiuned 
local aerodynamic centers based on calculations for Ordinary swept and 
imswept wings by means of the method of reference 2. 


The aerodynamic centers calculated in both ways have been expressed 
as fractions of the mean aerodynamic chord rearward of the leading edge 
of the mean aerodynamic chord. For wings with linearly varying chords, 
that is 

c = cj,[l - (1 - X)y^ (7) 


the mean aerodynamic chord is equal in size to the chord at the centroid 
of wing area defined by . . .. 


V* = i 

3 1 + ?^ 



and, hence, is equal to 


_ 2 i + X + 

- ""r 3 — iTx 



The spanwise location of this chord is either in the; plane of symmetry, 
when the wing as a whole is considered, or at the geometric centroid of 
wing area y, if one wing only is considered. The longitudinal locatlon_ 
of the mean aerodynamic chord coincides with the chord ^ the centroid 
of wing area in the case of unswept and ordinary swept wings. In the 
case of M, W, and A wings the quarter-chord point of the_mean aerodynamic 
chord is located at the distance ~ 


d = y tan 


tan Aq 
tan Aj^ 



^ y* \ 

1 + 2X 


( 10 ) 


rearward of the Intersection of the quarter-chord line and the plane of 
symmetry, where Aq and Aj_ are both positive for swee^ack. 


(W.BWT.IAIi^ 
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Scope of Calculations 

The plan foims for which calculations have been made are listed in 
table r as wings 1 to 20; they include eight M wings, eight W wings, 
three A wings and one inverted A wing. (The suggestive designation 
V wing is not used for the inverted A wing to avoid confusion with the • 
usage of that teim in connection with ordinary swept wings. ) All wings 
have a taper ratio of l/2; all have angles of sweep of either zero 
or except wings 10 to 15, which have angles of sweep of ±30°; all 

wings have an aspect .ratio of 6 except wings l6 to 20, which have an 
aspect ratio of 12. Three values of the spanwlse position of disconti- 
nuity in sweep, hereafter referred to as the "break, " are Included in 
this series of plan forms, namely y*g = 0.3^ 0.5^ 0.7* 

’For all plan foims lift distributions were calculated for unit angle 
of attack across the span, for linear symmetrical twist from 0 at the 
root to unit angle of attack at the tip, for" linear antisymmetric twist 
with unit angle of attack at the tip, and for unit effective angle of 
attack due to deflection of a 50-percent- span outboard aileron (ag& = l) . 
These lift distributions were Integrated to obtain total lifts, rolling 
moments, induced drags, and wing aerodynamic centers. 

For the sake of comparison lift distributions for a unit angle of 
attack across the span and for unit linear antisymmetric twist have also 
been calculated for eight ordinary swept and unswept wings, which are 
listed in table X as plan forms 21 to 28. The method of reference 1 was 
used in these calculations. Rolling-moment coefficients due to aileron 
deflection have not been calculated for the ordinary swept and unswept 
wings; they may be obtained from reference 5- However, only three 
points on the semispan and three terms in the series development for 
the lift distribution were used in the calculations of reference 5 ■ 
Therefore, the results of these calculations may not be as accurate 
as those of this paper for wings with large values of the parameter 
A/cos a. 


RESULTS AND DISCUSSION 
Spanwlse Lift Distribution 


The calculated spanwlse lift distributions for the twenty M, U, 
and A plan forms are presented for all four unit angle -of -attack condi- 
tions considered in this paper in figures 3 to 22. The symmetrical lift 
distributions are plotted as cc ^ JcC j,^ against the dimensionless span- 

wise ordinate y*; the antisymmetric lift distributions, as ccj/cC^ 
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against y*. In order to obtain values 'of the loading coefficient ccj^'C 

for the four angle-of-attack distributions these functions must be multi- 
plied by the values of and by the value of ^ "'*^Lp) P^sented in 

table T, respectively. — 

In order to afford a comparison of the lift distributions of these 
wings with each other and with lift distributions of ordinary swept and 
unswept wings, the spanwise lift distributions for constant unit angle 
of attack across the span are presented also in figures. 23 to 25* 

Effect of sweep on spanwise lift distribution .- The lift distri- 
butions of four M and W wings of aspect* ratio 6 with y* = 0.5 are shoTm 
in figure 23(a) with the angle of sweep of the inner portion of the wing 
as a parameter. ■ Also shown are the lift distributions of four ordina:^ 
swept wings of aspect ratio 6 with angle of sweep as a parameter. The 
curve labeled = -U5° pertains to plan fom 2, an M_jwing the inner 
and outer parts of which are swept -45° and respectively; similarly, 

the curve labeled = -30° refers to plan form 11, an M wing swept 
±30°, and the ones labeled A^ = 3^° -^i = ^5° 1 re?or to plan ^ 

forms l4 and 5., which are W wings swept ±30 and ±^5°^ respectively. The 
curve labeled A^ = 0 pertains to an unswept^wlng of ^pect ratio 6, 

which may be considered to be the limiting case of an M..pr a W wing as 

the angle of sweep of both the inner and outer parts of._the wing 
approaches 0. ^ 

The lift distributions of the ordinary swept wings shown in fig- 
ure 23 (^) exhibit the well-known characteristics of ’ sucE wings. As the”' 
angle of sweep increases positively from a sweptforward to a sweptback ^ 
wing the peak of the lift distribution at the root diminishes and turns 
into a dip; at the same time the lift near the tip increases. A similar^ 

behavior is noted for the M and W wings. The lift distrtbution at the 

roots of the wings exhibits a peak when the inner portion is swept^ for-_ 

ward and a dip when it is swept back; however, neither is as pronounced 

as the ones of the ordinary swept wings. The region near the breafc of 
the M and W also acts in a manner similar to that of the root of an 

ordinary swept wing; when the two halves of the wing ar^ swept forward 

relative to the break, as in the case of a W Wing, the ^ft distribution 
peaks near the break, whereas in the case of an M wing l,t dips at the _ 

break. Outboard of the break the M and ¥ wings act ;as if they were 
sweptback and sweptforward wings, respectively, with their roots at" the 
break, as may be expected, with a relatively rapid drop, to zero in the _ 
case of the W wings and a relatively high level of ilft^e.ar the tip in 
the case of the M wings. 

A set of lift distributions similar to those shown._for wings of 
aspect ratio 6 in figure 23(a) is shown in figure 23(b)^for wings of 
aspect ratio 12; however, lift distributions are shown dnly for M and W _ 
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and ordinary swept wings with or A = ±k^° and for an ordinary 

unswept wing. The curves are very similar to those of figure 23 (a) ^ 
except that for the wings with the higher aspect ratio the peaks are 
more pronounced and, in some cases, somewhat more localized than they 
are for those with the lower aspect ratio. 

Effect of spanwise location of break on spanwlse lift distribution . - 
The effect of the position of the discontinuity in sweep angle on the 
spanwlse lift distrihutlon is shown in figure 2^. The lift distributions 
for plan forms 1, 2, 3 (M wings of aspect ratio 6 with sweep of ±45°) as 
well as of plan forms 21 and 25 (ordinary swept wings of aspect ratio 6 
with 45° sweepforward and sweepback, respectively) are shown in fig- 
ure 24(a). The ordinary swept wings are Included because they constitute 
the limiting cases of M wings as the position of the break approaches the 
tip and the root, respectively. Figure 24(a) indicates that as the posi- 
tion of the break moves from the root to the tip the dip associated with 
the break moves outboard with the break; .also, a peak appears in the lift 
distribution at the wing root and becomes higher and less localized as 
the position of the break approaches the tip. Similarly, the lift distri- 
butions of ¥ wings of aspect ratio 6 and ±45° sweepback, plan forms 4, 5^ 
and 6, are shown in figure 24(b) along with the same ordinary swept wings 
as in figure 24(a). Figure 24(b) indicates that as the position of the 
break moves outboard the dip in the lift distribution becomes lower and 
generally less localized and that the peak associated with the break 
moves outboard with the break. 

Figures 24(c) and 24(d) are similar to figures 24(a) and 24(b), 
except that they represent M and ¥ wings of ±30°. The limiting cases as 
the position of the break approaches the root or tip are plan forms 22 
and 24. Although the peaks and dips of the lift distributions of fig- 
ures 24(c) and 24(d) are less than those of figures 24(a) and 24(b) as a 
result of the smaller angles of sweep, the peaks and dips associated with 
the break and their movement outboard with the portion of the break are 
just as distinct. Figures 24(e) and 24(f) are also similar to flg- 
xires 24(a) and 24(b), except that they represent wings of aspect ratio 12 
so that the limiting cases as the position of the break approaches the 
root or tip are plan forms 26 and 28. Again, although the peaks and dips 
are more pronomced than those observed in figures 24(a) and 24(b), the 
general trends of the lift distributions are the same. 

Effect of sween of outer nortion on snanwise lift distribution . - 
The effect on the lift distribution of the angle of sweep of the outer 
portion of the wing at a given angle of sweep of the inner portion of 
the wing is shown in figure 25. The three plan forms represented in 
figure 25(a) are wing 25^ an ordinary sweptback wing, wing 8, a A wing, 
and wing 6, a ¥ wing; all three plan forms have an aspect ratio of 6 and 
are swept back 45° iu the inner JO percent of the semi span. A comparison 
of the three lift distributions indicates that as the angle of sweepback 
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of the outer 30 percent of the wing is decreased the lift increases 
slightly near the tip and decreases slightly near in the inner portion ~ 
of the wing, with a peak near the break. As the outer ^portion is swept 
forward the peak beccones more pronounced, the lift near** the tip falls 
off, and the lift near the root is increased somewhat. A series of plan 
foims similar to those of figure 25(a) but with the: position of the bre^ 
at 30 percent of the semi span is represented in figure 25(b). For the _ 
plan forms represented in figure 25(b) the peak in the Tift distribution 
becomes more pronounced, the lift at and near the rOot Increases, and 
the lift near the tip decreases steadily as the outer 70 percent of the 
span is swept forward from a sweptback to a sweptforward position. 

The series of three plan forms represented in flgi^e 25(c) consists 
of an M, and inverted A, and an ordinary sweptforward wing of aspect 
ratio 6 with the inner 70 percent of the semispan syept forward ^5°* 
Comparison of the three lift distributions indicates that as the an^e 
of sweepforward of the outer 30 percent of the wing; is decreased and^ ~ " 
then changed to sweepback the lift near the tip increase steadily, but 
at the break it increases at first and then decreases ^i^tly, wEereas 
et the root it decreases sll^tly at first and then' increases quite 
rapidly. 

Figure 25 (<i) pertains to a series of plan forms l^ntical to that 
of figure 25(a), except that the wings have ah aspebt ratio of 12 rather 
than 6. The trends discussed in connection with fi^re"' 25(a) are more 
pronounced in the case of the series of plan forms represented in fig- 
ure 25 (d), but otherwise the same. ” 

Summary of effects of various narameters on snanwise lift distri- 
bution . - As a result of these considerations of thej 11§ distributions 
of M, W, and A wings in comparison with those of ordinary swept and 
unswept wings the conclusion may be drawn that the lif-^distributiohs ~ 
of the M, W, and A plan forms exhibit the characterjlst0s which may'Ee“ 
expected qualitatively if they are considered to be' the composites of 
ordinary wings; for instance, the break and the outer portion of an M 
wing tend to act as the root and one wing of a sweptback wing, whereas 
the break and the outer portion of a W wing tend to exhibit the charac- 
teristics of an ordinary sweptforward wing. The effec^ of the bre^ 
are more pronounced for wings with larger angles of' sw^p than for those 

with smaller ones. These effects also tend to be mbre 'pronounced, 

although sometimes more localized, in the case of wings with relatively 
high aspect ratios than for wings with moderate or lowPTspect ratios', as 
may be expected from a knowledge of the lift distribution of ordinary 
swept wings. 
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Lifts, Moments, Induced Drags, and 
Aerodynamic Centers 

The spanwlse lift distributions discussed in the preceding section 
have been integrated to obtain lift-curve slopes, coefficients of damping 
in roll, spanwise centers of pressure of the additional lift distri- 
butions, induced-drag coefficients, rolling-moment coefficients due to 
unit equivalent aileron deflection and aerodynamic centers. The 

aerodynamic centers have been computed both by assuming the local aero- 
dynamic centers to be at the quarter-chord line and the somewhat arbi- 
trarily chosen locations shown in figure 26, which are based on the 
results of calculations by the method of reference 2. The results are 
presented in table I and figures 27 to 29 for the M, W, and A wings con- 
sidered in this paper as well as for the comparable ordinary swept and 
unswept wings. The wing aerodynamic centers obtained by means of the 
local aerodynamic centers given in figure 26 will be referred to as the 
corrected aerodynamic centers and are listed as (a. in table I 
and figures 27 to 29. 

Effect of sweep on lift, drag, and aerodynamic center . - The 11ft- 
curv© slopes, induced-drag ratios and the aerodynamic centers 

of M and W wings are shown plotted as fxmctlons of the angle of sweep of 
the inner wing portion with the position of the break as a parameter in 
figures 27(a) and 27(11) for wings of aspect ratios 6 and 12, respectively. 
The curves in figures 27(a) and 27(b) for y*^ = 0.5 pertain to the 

series of plan forms represented in figures 23(a) and 23 (t), respec- 
tively. The curves for y*_ = 0 and y*_ = 1 represent ordinary 

Jd J5 

swept wings which constitute the limits of M and W wings as the break 

approaches the root or the tip. Consequently, whereas the curves 

for y*2 = 0.3, 0.5^ and 0.7 represent wings which change frcm an M wing 

to an ordinary unswept to a ¥ wing as the angle of sweep of the inner 

portion of the wing increases from -i^5° through 0° to the curve 

for y*g = 0 represents an ordinary wing which varies from the 45° swept- 

back through the unswept to the 45° sweptforward position, and that 

for y*T3 = 1, an ordinary swept wing which changes from the 45° swept- 
a 

forward to the 45° sweptback position. 

The lift-curve slopes of all plan forms represented in figure 27 
are highest for the unswept position and decrease as the plan forms are 
swept either backward or forward, either into M and ¥ or into ordinary 
swept wings. The decrease in lift-curve slope is less for the true M 
and ¥ wings with the breaks at 0.3, 0.5^ and 0.7 of the semispan than 
for the ordinary swept wings or M and ¥ wings with the breaks very near 
the root or tip. The Induced drag of the M wings is slightly higher than 
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that of the ordinary swept wings , but that of the W wings with the 
break at 0.3 and 0.5 of the semispan is slightly lower than that of the 
ordinary swept wings for angles of sweep between about 10° and 35°* The 
aerodynamic centers are more rearward for both the swept and the M and 
W wings than for the unswept wing. The aerodynamic centers of the W wii^s 
are between the limits of the equivalent sweptback and sweptforward wings, 
but the aerodynamic centers of the M wings are farther forward than those 
of the equivalent ordinary swept wings. These trends are true both for 
the uncorrected and the corrected aerodynamic centers. '".'Both are in goo3. 
agreement with each other, within 1 or 2 percent of .chord for the M and_ 

W wings and within 4 percent for the ordinary swept wings; the corrected 
values are generally slightly farther forward than the uncorrected ones. 

For the wings of aspect ratio 12 represented in figure 27(b) the 
lift-curve slope decreases with sweep in almost the same manner for M 
and W as for ordinary swept wings. The induced drag of M wings is 
larger than that of the ordinary swept wing, whereas that of W wings is 
the same as that of the ordinary swept wing. The uhco:^^cted aerodynamic 
centers (no corrected ones having been calculated) move rearward as the 
plan forms are swept, the change being less for the W wings than the 
ordinary swept wings and still leas for M wings. Tlie aerodynamic centers 
for y *2 = 0.3, 0. 5.> aud 0.7 are generally within about 0.02 chord of 

each other. ■■ - ' f „'. 

Effect of snanwise location of break on the lift, -drag, and aero- 
dynamic center .- Figures 28(a) and 28(b) constitute! cnJgs plots of fig- 
ures 27(a) and 27(b), respectively. The aerodynamifc p^ameters are 
plotted against the position of the break with the angK of sweep as a 
parameter. In effect, each of the series of plan forms; for which lift 
distributions are shown in figure 2k is represented; by a curve (for each 
of the aerodynamic parameters) in figure 28. The curves ' of figure 28 ~ 

serve to corroborate the conclusion reached in exam^nii^ figure 27^^^ 
the effect of the position of the break is relatively small on both "the 
lift-curve slope and on the induced-drag coefficient; fn the case of the 
wings with taper ratio 0.5 considered in this paper, the" least induced 
drag is Incurred for W wings with angle of sweep betwe^ 0° and 30° and 
the position of the break between 0.2 and 0,6 of th^ s^iispan. Except 
for the M wing with i<-5° of sweep the aerodynamic cehter" moves generally 
forward as the position of the break is moved outboard from the root to" 
the tip. 

Effect of sween of outer TJortion on the lift, drag, and aero- 
dynamic center .- The effect on the aerodynamic coefficients of changes 
in the angle of sweep of the outer wing portion while maintaining that 
of the inner wing portion at a given value is shown injfigure 29, which 
represents series of plan forms similar to tHose represented in fig- 
ure 25 . If the position of the break, which is the paFameter of the 
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curves of figure 29 , is at the root, the entire ving is swept in these 
series; if the position is at the tip, the plan form is unchanged and is, 
in the case of figures 29 ( b .) and 29(c), a svepthack wing and in the 
case of figures 29(13) and 29(d), a 45° sweptforward wing. If the posi- 
tion of the break is at 0. 3, 0. 5^ or 0. 7 of the semi span the plan form, 
in the case of figures 29(a.) and 29(c), changes from a W to a A to an 
ordinary sweptback wing as the angle of sweepback of the outer panel is 
increased from -45° through 0° to 45°; in the case of figures 29(b) and 
29(d) the plan form changes from an ordinary sweptforward to an inverted 
A to an M wing as the angle of sweepback is increased from to 45°. 

Figure 29(a) shows that as the angle of sweepback of the outer por- 
tion is increased from -45° to 45° the lift-curve slope first increases 
then decreases, with a maximum near the true A-wing condition of an 
unswept outer wing portion. This maximum is higher when more of the 
wing portion is involved in this sweeping process, that is, the closer 
the position of the break is to the wing root. The induced drag is a 
minimum for the ordinary swept wing represented by y*^ = 0 between 0° 

and 10° of sweepforward. For y*^ =0.5 and 0.7 the Induced drag is 

hipest near the A-wing condition and less for either the W-wing or 
ordinary sweptback-wing condition; for y*^ * 0. 3 the drag appears to 

be substantially independent of the angle of sweep of the outer panel. 

The aerodynamic centers of the plan forms with y*^ = 0.3 and 0. 5 are 

between those with = 0 and = 1.0; the curve for y*^ = 0.7; 
however, is on the side of the curve for y*^ =1.0 away from that 
for y*g = 0. For values of y*^ = 0, 0.3, and 0.5 the aerodynamic 
center is most forward for the A-wing condition, for y*^ =0.7 It is 
most rearward for this condition, and for y*^ =1.0 it is constant, 

since the plan form is the same for all values of the ordinate, Aq. As 
noted in connection with figure 27(a) the corrected and izncorrected 
aerodynamic-center values are generally within 0.02 chord of each other, 
the corrected ones being slightly more forward on the average. 

For the wings represented in figure 29(b) with the inner portion 
swept forward 45°, the lift-curve slope has a metximxam when the outer 
portion is unswept (the Inverted A-wing condition) , the maxi m u m value 
being higher the farther the break is Inboard. ‘ The induced-drag curves 
have minima between 0° and 15° sweepforward of the outer portion, the 
drag being the lower the more of the wing is swept back, that is, the 
farther the break is inboard. The aerodynamic centers move forward as 
the outer panel is swept back from -45° to 45° of sweepback, except that 
if the entire wing is swept, that is, if y*^ = 0, the aerodynamic center 

moves rearward for positive values of Aq. The corrected aerodynamic 
centers are generally within about 0.03 chord of the uncorrected ones 
and slightly farther forward on the average. 
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Figures 29(c) and 29(d) are similar to figures . 29(a) and 29(13) but 
pertain to vings of aspect ratio 12. The trends of curve 8 0? the 
aerodynamic coefficients as functions of the angle Of sweepback of the 
outer wing portion are very similar to those shown in figures 29(a) 
and 29 (b ) . 


Limitations of the Results 

The lift distributions calculated in this paper as well as the 
other aerodynamic parameters obtained by integrating the lift distri- 
butions are subject to certain limitations in consequence of the approxi- 
mations made in the calculations. These approximations may be considered 
to be of two types, those in the manner in which a true lifting-surface 
(thin-wing, potential-flow) solution is effected by th^ method of calcu- 
lation used in this paper, and those in the degree to which a thln-wing' 
potential-flow solution represents reality. 

The quarter-chord-vortex three-quarter-chord-downwash concept is 
known to be capable of furnishing generally excellent results for the 
spanwise lift distribution, although it furnishes no information what- 
ever concerning the local centers of pressure of lift distribution. 

This concept fails near the root of swept wings, since'^n this region 
the chordwise pressure distribution cannot be approixlm^ed by a super- 
position of two-dimensional angle-of-attack-type and parabolic-camber- 
type pressure distributions. For ordinary swept wings "Ihis shortcoming*’ 
is not serious, and its effects, which appeaF to be smallj, can be local- 
ized by using many points along the span in methods' of lihe type of those 
of references 1 and 2 or by considering each vortex individually as in 
the present paper. For an M, ¥, or A wing, however, the effects of this 
shortcoming of the one -quarter -chord-vortex three -quar'fer-chord-downwash 
concept are somewhat more serious, because it affects ffie regions near ' 
the break as well as near the root. 

Another limitation on the results calculated in this paper is 
imposed by the deviation of potential theory from actuality. On swept ' 
wings the boundary layer tends to flow along the span and accumulate in 
the downstream, regions of the wing, that is, the tip o^a sweptback wing 
and the root of a sweptforward wing. This accumulation causes a decrease 
in the lift in the affected region. Consequently, pweptforward wings do 
not always exhibit the peak in the lift distribution predicted by theory, 
nor do sweptback wings carry as high a level of lift near their tips. 

The same phenomenon must be expected on M and W wings. Consequently, 
the magnitude of the peak in the lift distribution at the break of a 
W wing may be reduced and, in some conditions, changed into a dip. 
Similarly, the M wing is likely to carry less lift near its ro'ot and ' 
tip than predicted by theory as a result of boundary- layer effects. 

These changes in the spanwise lift distribution are likely to be reflected 
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in a lower lift-curve slope, a more forward aerodynamic center than 
given in this paper, and an induced drag which may be either somewhat 
larger or somewhat smaller than that predicted on the basis of the calcu- 
lations of this paper, depending on the nature of the change in the lift 
distribution. 

Inasmuch as the theoretically computed results are qualitatively 
in line with the ones that might be expected from experience with ordl- 
naiy swept wings, there is a possibility that a knowledge of the effects 
of the boundary layer on the aerodynamic characteristics of ordinary 
swept wings may seive to correct the results presented in this paper for 
such effects. Another possibility is that the accumulation of the 
boundary layer may be controlled, by means of suitable fences, for 
Instance, in which case the results presented in this paper would be 
more nearly applicable. 

The calculation of the aerodynamic centers in this paper is open to 
question on the grounds that the local aerodynamic centers are not known, 
those presented in figure 26 being in the nature of carefully considered 
estimates only. However, comparison of the wing aerodynamic centers 
calculated on the basis of these values with those calculated on the 
assumption of local aerodynamic centers at the quarter chord are in good 
agreement, those calculated by means of the estimated aerodynamic centers 
being only slightly farther forward on the average. Consequently, the 
aerodynamic centers are probably as accurate as the other aerodynamic 
parameters presented in this paper. 

Compressibility effects have not been taken into account in the 
calculations of the lift distributions and the other aerodynamic param- 
eters. With angles of sweep between 30° and ccmpressibility effects 
are not likely to be Important up to free-stream Mach numbers of about 
0.6 or 0.7. At higher free-stream Mach mmubers the three-dimensional 
Glauert-Prandtl correction may serve to correct the results presented 
in this paper for compressibility effects in an approximate manner. At 
Mach numbers higher than 0.8 or O.p, however, the corrected results are 
likely to possess only limited utility. 


CONCLUDING RMAEKS 


The results of theoretical incompressible-flow calculations of the 
spanwlse lift distributions, lift-curve slopes, spanwlse centers of 
pressure, aerodynamic center, coefficient of damping in roll, and 
rolling-moment coefficient due to aileron deflection of twenty M, W, 
and A plan forms have been presented and compared with similarly calcu- 
lated results for ordinary swept and unswept wings. The aerodynamic 
characteristics of M, N, and A wings are, qualitatively, those that 
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■would be expected on the basis of knowledge of the characteristics of 
ordinaiy swept wings. The theoretically calculated induced drags of 
these wings differ little from those of ordinary swept wings. 


Langley Aeronautical Laboratory 

National Ad'vlsory Committee for Aeronautics 
Langley Field, Va. 
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TABLE I.- BUMMAKY OF OKMETHIC AHB ASiOIlYHAMIC PARAHETERfi - Concluded 



1 Qe.ometric parameters 

Aerod^amlc parameters 

Wing 



'' 













Plan fora 

Type 

A 

\ 

(ieg) 

(deg) 




-C7 

y* 

a. c. 

(a. c. 


15 


w 

6 

1/2 

30 

-30 

.7 

4.13 

-0543 

.421 

.444 

.260 


.459 

l6 


M 

12 

1/2 

-45 

45 

.3 

3.91 

.0350 

.472 

.450 

.322 


.480 

17 


K 

12 

1/2 

-45 

45 

.5 

3.85 

.035? 

.477 

.438 

.326 


.537 

18 

w 

W 

12 

1/2 

45 

-45 

.5 

4.02 

.0294 

.491 

.448 

.408 


.536 

19 


W 

12 

1/2 

45 

-45 

.7 

3.87 

.0295 

.483 

.k6k 

.4U 

• 

.521 

20 


A 

12 

1/2 

45 

0 

.7 

4.07 

.0318 

.530 

483 

.277 


.578 

21 

V 

Svept 

6 

1/2 

-45 

-45 


3.33 

.0575 

.342 

.378 

,418 

.371* 


22 


Swept 

6 

1/2 

-30 

-30 


3.91 

.0546 

.388 

.401 

.318 



23 

C 1 

Unswept 

6 

1/2 

0 

0 


4.32 

.0530 

.413 

.425 

.250 

.239 


2h 


Swept 


1/2 

30 

30 


3.98 

.0540 

.390 

.444 

.255 



25 


Swept 

6 

1/2 

45 

45 


3.44 

.0556 

.355 

.455 

. 29 ^ 

.310 


26 


Swept 

12 

1/2 

-45 

-45 


3.87 

.0304 

.469 

.341 

.700 



27 


Unswept 

12 

1/2 

0 

0 


5.16 

.0276 

.555 

.429 

.250 



28 


Swept 

12 

1/2 

45 

45 


3.94 

.0299 

.457 

MS 

.470 






f 



Corrector 

'vortex 
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Downwa-^h point 


Figure 1.- Location of horseshoe vortices aad.dovnwash points on an 

arhltrary plan form. ' 
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(a) Symmetrical lift distributioiis. 

(b ) Antlsymmetrlcal lilt distributions. 

Figure 3.- Spanwlse lift die trlbut ions for plan form 1 (A = 6, A, = 0.5, 

A = M wing, jtb* = 0.3). 


DimenslonleBS spanwlse ordinate, /■ 

(a) Symmetrical lift distributions. 

^(b) Antlsymmetrical lilt distributions. “ 

Figure 4.- Spanwlse lift distributions for plan fom ^ (A = 6 

A= M wing, y^* = 0.5); “ 







Spanwlse loading coefficient, ^ Spanwlse loading coefficient 
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Dimensionless spanwise ordinate, y* 

(a) Symmetrical lift distributions. “ 

(b) Antisymmetrical lUt distributions. 

Figure 8.- Spanwise lift distributions for plan form 6" (A = 6, X 

A = ±45°, W wing, y^* =0.7). 
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(a) Symmetrical lift distributions, 
(b) Antisymmetrical lift distributions. 


Figure 9 .- Spartwlse llffc distributions for plan fom J (A = 6, \ = 0. 5, 

\ = ^5°, Awing, yB* = 0.3). 
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(a) Symmetrical lift distributions, 
(b) Antisymmetrical lift distributions. 


Figure 11.- Spanwlse lift distributions for plan form 9 (A 
X = 0.5j ^ Inverted A ving, = O.T). 


.8 .9 1.0 

Dimensionless spanwise ordinate, y* 

(a) Symmetrical lift distributions, 

(b) Antlsymmetrical lift distributions. 

Figure 12.- Spanwise lift distributions for plan form 10 (A = 6, 

X = 0.5, A = ±30°, M wing, y^* - O. 3 ). 
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Dimensionless spanwise ordinate, y* 

(a) Symmetrical lilt distributions, 
(b) Antisymmetrical lift distributions. 


Figure 13.- Spanwise lift distributions for plan form 11 (A = 6 
A, = 0.5, A = ±30°, M wing, = 0.5). 
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(a) Symmetrical lilt distributions, 
(b) Antlsymmertrlcal lift distributions. 


Figure 15.- 


Spanwise lift distributions for plan form 13 (A 
1 = 0.5, A = ±30°, W wing, 7 b* = 0.3). 
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Spanwise loading coefficient, - Spanwise loading coefficient, r 
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(a) Symmetrical lift distributions, 
(b) AntlEyiumetrloal lift distributions. 


Figure 21. - Spanwlse 
X = 0.5, 


lift distributions for plan foha I9 (A - 12 
A= ±11.50^ w wing, yg* = 0 . 7 ). 
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Figiire 22.- 


Spanvlse lift distributions for plaji foinn 20 (A = 12, 
X = 0.5, Ai = 45°, A wing, 73 * = 



ra 


T 


.5 


.6 


Dimensionless spcmise ordinate^y* 


(a)/\=6. 

Figure 23.- Effect of angle of sweep on the lift dlatrlbutlona 
W, and ordinary swept wings of taper ratio 0.5. 
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{f)WmngSi A=i2,yip4S°. 


Figure 2k.- Concluded. 
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Anqle of sweep of the inner portion of the winq,wi; , deg 

(b)A=l2. ■ y - 

Figure 27. - Concluded. 
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Dimensionless ordinate of point of dlscontinulty^yQ 

Figure 28.— Effect of location of the spanvlse discontinuity In the 
angle of sweep on the aerodynamic characteristics of M and W wings, 
taper ratio 0.5. 
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Figure 29.- Effect of sweep of outer portion on aerodynamic 

characteristics , 
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Angie of sweep of theater portion of the wing^wi^^deg 


Figure 29 .- Continued, 
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Anqle of sweep of the outer portion of the winq,uirr>deQ 

(d)yi/.-45‘’,A^l2. ' r ^ 

Figure 29.- Concluded. 
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